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ABSTRACT 

Ultra-high cosmic rays (UHECRs) with energies > 10 19 eV emitted at cosmological distances will be 
attenuated by cosmic microwave and infrared background radiation through photohadronic processes. 
Lower energy extra-galactic cosmic rays (~ 10 18 — 10 19 eV) can only travel a linear distance smaller 
than ^Gpc in a Hubble time due to the diffusion if the extra-galactic magnetic fields are as strong as 
nano Gauss. These prevent us from directly observing most of the UHECRs in the universe, and thus 
the observed UHECR intensity reflects only the emissivity in the nearby universe within hundreds of 
Mpc. However, UHECRs in the distant universe, through interactions with the cosmic background 
photons, produce UHE electrons and gamma-rays that in turn initiate electromagnetic cascades on 
cosmic background photons. This secondary cascade radiation forms part of the extragalactic diffuse 
GeV-TeV gamma-ray radiation and, unlike the original UHECRs, is observable. Motivated by new 
measurements of extragalactic diffuse gamma-ray background radiation by Fermi/LAT, we obtained 
upper limit placed on the UHECR emissivity in the distant universe by requiring that the cascade 
radiation they produce not exceed the observed levels. By comparison with the gamma-ray emissivity 
of candidate UHECR sources (such as GRBs and AGNs) at high-redshifts, we find that the obtained 
upper limit for a flat proton spectrum is ~ 10 1 ' 5 times larger than the gamma-ray emissivity in GRBs 
and ~ 10 times smaller than the gamma-ray emissivity in BL Lac objects. In the case of iron nuclei 
composition, the derived upper limit of UHECR emissivity is a factor of 3-5 times higher. Robust 
upper limit on the cosmogenic neutrino flux is further obtained, which is marginally reachable by the 
Icecube detector and the next-generation detector JEM-EUSO. 
Subject headings: cosmic rays- gamma-rays: diffuse background-neutrinos 
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1. INTRODUCTION 

There is a general consensus that cosmic rays with 
energy above 10 19 eV originate from extragalactic 
astrophysical sources, although the sources are uniden- 
tified. Some candidates have been proposed, including 
active galactic nucleus (AGN) jets (e.g. Biermann & 
Strittmattcr 1987; Bcrczinsky ct al. 2006), gamma-ray 
bursts (GRBs) (Waxman 1995; Vietri 1995; Wick et 
al. 2004; Murase et al. 2008), and semi-relativistic 
hypernovae remnants (Wang et al. 2007). Any viable 
candidates must be able to provide a right amount of 
UHECR emissivity to match the observed flux. Due 
to the attenuation by cosmic microwave and infrared 
background radiation through photohadronic processes, 
the energy loss distance of UHECRs above 10 19 eV is 
less than several hundreds of Mpc. Although extra- 
galactic cosmic ray protons in the lower energy range 
10 18 < E < 10 19 eV suffer from much less energy 
loss, they could possibly travel a linear distance much 
smaller than the size of the universe as well due to 
the diffusion in the intergalactic magnetic fields (e.g. 
Lemoine 2005). For UHECRs whose Larmor radius Vl 
is much larger than the coherence scale (l c ) of the fields, 
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which is valid when E > 10 18 eV(£/lnG)(Z c /lMpc), 
the scattering length is i S catt — r \l^c and thus 
the diffusion coefficient^ is D{E) — (l/3)cZ sca tt = 
4 x 10 34 ( J B/10 18 eV) 2 (S/lnG)- 2 (Z c /lMpc)cm 2 s- 1 . 
In a Hubble time, these particles travel 
a linear distance d ~ (2Hq 1 D) 1 /' 2 ~ 
60( J B/10 18 eV)( J B/lnG)- 1 (Z c /lMpc)~ 1 / 2 Mpc, where 
Hq is the Hubble constant at present. Therefore, the 
observed UHECR intensity reflects only the emissivity 
in the nearby universe within several hundreds of Mpc. 
Usually, in the literature this emissivity is compared 
with the emissivity in electromagnetic radiation pro- 
duced by local candidate sources to check whether 
they can provide sufficient power. The same type 
of candidate sources are also present in the distant 
universe and produce electromagnetic radiation and 
UHECRs as well. It would be also useful to compare the 
emissivity of UHECRs and emissivity in electromagnetic 
radiation in the distant universe, especially when the 
local sources of the candidates, such as GRBs or AGN 
giant flares (Farrar & Gruzinov 2009), are too infrequent 
to be detectable in years. While the emissivity in the 
electromagnetic radiation can be measured directly for 
distant sources, UHECRs produced by them can not 
be observable and therefore their emissivity can not 
be determined directly However, UHECRs produce 



5 For extragalactic cosmic ray protons with E < 
10 18 eV(B/lnG) _1 (Z c /lMpc), the smaller angle diffusion ap- 
proximation will not be valid and the diffusion coefficient in this 
regime will be different. 
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UHE electrons and gamma-rays that in turn initiate 
electromagnetic cascades on extragalactic background 
light (EBL) and this lower-energy, secondary cascade 
radiation can be observable. Thus the energy of UHE- 
CRs will be eventually converted to observable diffuse 
gamma-rays, as was first noticed by Wdowczyk et al. 
(f 972). The spectrum of this cascade radiation is rather 
insensitive to the spectrum of the original UHECRs 
(e.g. Strong et al. 1973; Berezinsky & Smirnov 1975; 
Coppi & Aharonian 1997), and thus the total level 
of the cascade background acts as a particle detector 
calorimeter, allowing us to measure the total UHECR 
cmissivity in the universe. 

Recently, Kalashev et al. (2009) studied the contri- 
bution of the cascade radiation by UHECRs interacting 
with EBL photons to the diffuse extragalactic gamma- 
ray background (EGBR) measured by EGRET. Using 
the new measurement of the EGBR by Fermi/LAT, 
Berezinsky et al. (2011) and Ahlers et al. (2010) at- 
tempted to constrain the UHECR source evolution mod- 
els and made predictions for the flux level of cosmogenic 
neutrinos. All these papers assume that the accumulated 
UHECR flux from a homogenous population of sources 
in the whole universe fit the observed energy spectrum of 
UHECRs. In our paper, in order to obtain an indepen- 
dent constraint on the emissivity of UHECRs in the dis- 
tant universe, we will relax this assumption since UHE- 
CRs above 10 19 eV are produced dominantly by local 
sources within hundreds of Mpc, much smaller than the 
size of the universe, whose density could well be enhanced 
or deficient relative to the average density. Indeed, there 
are suggestions that one single source in the nearby uni- 
verse, such as Cen A, is the dominated source producing 
the observed UHECRs (e.g. Cavallo 1978; Farrar & Pi- 
ran 2000; Rieger & Aharonian 2009). 

In §2, we first present the constraints on the cascade 
emission imposed by the Fermi/LAT observations of the 
EGBR. Then in §3.1, we present an analytic approach 
for calculating the energy density of the cascade emis- 
sion produced by UHECRs in the universe. By compar- 
ison with the allowed maximum energy density of the 
cascade radiation measured by Fermi/LAT, we obtain 
upper limits on the UHECR emissivity at high redshifts 
in §3.3. We then discuss the effect of synchrotron loss of 
secondary electrons in the presence of intcrgalactic mag- 
netic fields on the cascade radiation in §3.4. Upper limits 
on the cosmogenic neutrino flux are further obtained in 
§3.5. We discuss the case of heavy nuclei composition of 
UHECRs in §4. Finally, we give a summary in §5. 

2. CASCADE RADIATION 

UHECRs interact with EBL photons through photo- 
pion process or Bethc-Hcitler process and produce UHE 
electrons and gamma-rays. These UHE electrons and 
gamma-rays interact via Compton and pair-production 
process with soft photons in the CMB and radio back- 
ground. This will lead to the development of an electro- 
magnetic cascade, in which the number of electrons and 
photons increase quickly. The cascade proceeds until the 
energy of photons drops below the pair creation thresh- 
old. This process reprocesses essentially all the energy 
of UHE electrons and gamma-rays produced by UHE- 
CRs into lower-energy photons below ~ TeV which we 
can detect. The cascade low-energy photon spectrum is 



a universal spectrum, described by 

dn ( e' 1 - 5 for e 7 < e b 

de 7 { e 7 for e b < e 7 < e cut 

with a steepening at e 7 > e cut , where e cut is the ab- 
sorption energy of a cascade photon scattering on EBL, 
e 6 = (e cu t/lTeV) 2 GeV, and a 7 ~ 1.8 - 2 typically (e.g. 
Berezinsky & Smirnov 1975; Coppi & Aharonian 1997). 
The peak of the energy distribution in this spectrum is 
at e cu t, which is estimated to be ~ lOOGeV (Berezinsky 
et al. 2011). By requiring that this theoretical cascade 
spectrum touches the lower end of the error bars of the 
Fermi/LAT data (Abdo et al. 2010), Berezinsky et al. 
(2011) obtained an upper limit on the cascade energy 
density, i.e. 

^7^^ar = 5.8xl0- 7 cVcm- 3 . (2) 

where w 7 is the cascade energy density produced by any 
UHECR source population in the universe. 

3. CONSTRAINING THE UHECR EMISSIVITY WITH THE 
DIFFUSE GAMMA-RAY EMISSION 

3.1. Energy loss of UHECRs into the electro-magnetic 
component 

Since the cascade emission is mainly in the GeV-TeV 
range and has an almost universal spectrum, it is nu- 
merically much more efficient to calculate the total en- 
ergy density w 7 injected into the cascade and compare 
this value with the limit w™f x imposed by Fermi/LAT. 
The cascade energy density is the accumulation of the 
electro-magnetic (EM) radiation produced by the source 
population over the whole universe. One can calculate it 
by summing the contributions by individual sources that 
are generated at different cosmological epochs, i.e. 

W 7=/ J L p {t)Pv,em{E-; iZ {t))dtdE-; (3) 

where 

dN 

Lp (t) = h(z)E;^ (4) 

is the emissivity of protons per energy decade at energy 
E p at some cosmological epoch t, /?o,em is the fraction 
of the proton energy deposited into the EM component 
(photons and pairs) that we observed at present, n(z) is 
the comoving-frame number density of protons injected 
per unit time at time t = t(z), E^ is the proton energy at 

the source, and is the energy spectrum of protons. 

The energy of injected protons evolves with time as 

dE 

-^=E p H(z) + b(E p ,t) (5) 
where H{z) is the Hubble constant at time t = t(z), and 

7 [o"BH (C7)/bH (£7) + ^7(^7)^7(^7)] 

(6) 

is the total energy loss rate of protons of energy E p due to 
photopair and photopion interactions with EBL photons 
at some cosmological time t(z), T is the proton Lorentz 
factor, cjbh and the cross section for photopair 



3 



and photopion production respectively, /bh and / P7 are 
the fractions of energy loss of protons due to photopair 
and photopion interactions in one collision, and n 7 (e, z) 
is the number density of EBL photons of energy e at 
rcdshift z. We use the cross section in Chodorowski et 
al. (1992) for the Bethe-Heitler process and use the full 
photopion production cross section from the pion pro- 
duction threshold up to high energies as described in 
Miickc ct al. (2000). The EBL includes CMB photons 
and infrared-to-optical background photons (Finke et al. 
2010). The number density of infrared-to-optical back- 
ground photons at high-rcdshifts are taken from the data 
set onlin^B- 

From the above equations, one can solve E p {t). Then 
one can obtain the energy loss of protons exclusively into 
the EM component per unit time at some cosmological 
time i(z), 



dE„ 



dt 

+R em a 



M = 

p 7 (£ 7 )/p 7 



E p i 
21"-' 



<te 7 [CT BH 

(e 7 )] x E n 



(£ 7 )/bh(£ 7 ) 



f 

Je- 



/2T 



de- 



(7) 



where R em — 0.6 is the fraction of proton energy that 
goes into the EM component in the photopion channel 
(the other ~ 0.4 goes into the neutrino production chan- 
nel) (Engel et al. 2001). 

Then the energy lost into the EM component by a pro- 
ton of energy E£ (at the source) during the whole period 
from the injection time i(z) to the present time is 



i( E p^) E p = / 
JO 



f ^ ^ dEp^em(z^) 



p,em \ 

~~dT 



-dt = 



dE p ^ em (z) dt 
dt 



dz. 



dz 

The corresponding energy that remains at the present 
epoch after taking into account the redshift energy loss 
is 



Po,em(Ep,t)Ep 



t(z) 



1 dEp ern [z^) 



l + z 



dt 



-dt. 



(9) 



The values of f3 ern (E p ,t) as a function of the energy 
of protons generated at different redshifts are shown in 
Fig.l. For a source at z > 0.5, more than a half of the 
proton energy is lost into the EM cascade for protons 
with energy above 10 185 eV. 

3.2. The UHECR source density evolution 

The evolution of the UHECR emissivity per energy 
decade with redshift can be parameterized as 

S(z) 

JpV -~ 11 s( z = iy 

where L p (z = 1) = h(z = l)E p ^t, n(z = 1) is the 

comoving-frame number density of protons injected per 
unit time at z = 1 and S(z) is the source density at 
rcdshift z. We take three cases for S(z) in the following 
calculation, i.e. S(z) follows the star formation history 
(SFR), gamma-ray burst (GRB) rate and active galactic 
nuclei (AGN) rate in the universe respectively. 
We take the form of SFR from Yiiskel et al. (2008) 



L p (z) = L p (z = 1) 



(10) 




http:/ /www. phy.ohiou.edu/ finke/EBL/index.html 



(11) 



Recent analysis of the GRB redshift distribution as 
detected by Swift reveals that the GRB rate is enhanced 
at high redshift relative to SFR (Le & Dermer 2007). 
This may arise from some mechanisms (Kistler et al. 
2008), such as a GRB preference for low-metallicity en- 
vironments (Stanek et al. 2006; Langer & Norman 
2006). Following Yuksel & Kistler (2007), we assume 
Sqrb(z) oc (1 + z) 1a S SF h, which gives 



Sgrb(z) 



oc 



(l + z) 4 ' 8 , 
(l + z) 1 - 1 , 
(l + z)" 2 - 1 



z < 1 
1 < z < 4 
z > 4 



(12) 



As shown by Yuksel & Kistler (2007), this source den- 
sity evolution function is consistent with the evolution 
function obtained in Le & Dermer (2007). 

AGNs may have a similarly strong evolution with red- 
shift, as found in Hasinger et al. (2005) for different 
luminosity AGNs. Following Hasinger et al. (2005) and 
Ahlers et al. (2009), we take the form of 



Sagn(z) oc 



(l + z) 5 - , 
constant, 

1Q (2.7-*)_ 



z < 1.7 
1.7 < z < 2.7 
z > 2.7 



(13) 



3.3. Constraints on the emissivity of UHECRs in the 
distant universe 

Once we know /3o je m and the source density evolution 
function S(z), we can calculate the cascade energy den- 
sity w 7 , i.e. 



Po,e.m{E S p ,z)L p (z = 1) 



S(z 



S(z) dt 

i)Tz dE ^ 

(14) 



energy of extragalactic cosmic rays, z max is the maximum 
redshift of extragalactic cosmic ray sources, and dz/dt = 
H (l + z)[0 M (l + z) 3 + Oa] 1 / 2 (with n = 0.3, n A = 0.7, 
and Ho = 70kms -1 Mpc~ 1 ). We assume a power-law 
spectrum for extragalactic cosmic rays between E PtW \ n 
and -Bp, max, i-e. dNp/dEp 1 oc (£^)~ 7s . As (3 em becomes 
negligible at E p < 10 17 eV, E p ^ min is taken to be 10 17 eV 
in the following calculations, unless otherwise specified. 

Requiring cj 7 < o;™| x , one can obtain the upper limit 
on the UHECR emissivity at z = 1 for different source 
density evolution scenarios. We define the total UHECR 
emissivity as the integral of L p [z = 1) over the energy 
range from E P}min to E p ^ max , i.e. 



P(z = 1) 



L p {z = l)dE s p . 



(15) 



The upper limits of P(z = 1) for different spectral in- 
dex 7 ff are shown in Fig. 2 (the upper panel) with a fixed 
maximum proton energy E pmax — 10 21 eV (but differ- 
ent minimum energies) for different source density evo- 
lution scenarios. The values are in the range of a few 
10 45 to a few 10 46 ergMpc _3 yr _1 . The upper limit for 
the SFH case is the highest because of the slowest red- 
shift evolution and smaller /3o,em at lower redshifts. As 
a comparison, we also show the gamma-ray emissivity of 
GRBs and BL Lac objects at redshift z = 1. Though 
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the local gamma-ray emissivity of GRBs inferred by dif- 
ferent groups are different by about one order of mag- 
nitude (Guetta et al. 2005; Le & Dermer 2007; Wan- 
derman & Piran 2010), the inferred emissivity in GRBs 
at redshift z = 1 has a much smaller uncertainty, be- 
ing 1.1 x 10 44 (Ai/10s)ergMpc- 3 yr _1 in Guetta et al. 
(20 05), 1.8 x 10 44 (At/10s)ergMpc~ 3 yr" 1 in Wanderman 

& Piran (2010), and 2.5 x 10 44 (At/10s)ergMp C - 3 yr" 1 
in Le & Dermer (2007), where At is the mean dura- 
tion of long GRB in the explosion frame. By compar- 
ison with the upper limits on the UHECR emissivity, 
we find that the ratio between the two emissivities is 
-PcR.(10 17 -10 21 eV)/P 7 < lO 15 for a flat cosmic ray spec- 
trum -fg ~ 2.0. 

Dermer & Razzaque (2010) recently obtained the 
gamma-ray emissivity of different classes of AGNs. The 
two bright classes at high redshift are BL Lac ob- 
jects and flat spectrum radio quasars (FSRQs). The 
gamma-ray emissivity in BL Lac objects and FSRQs 
at z = 1 are, respectively, 6 x 10 46 ergMpc _3 yr _1 and 
1.5 x 10 46 ergMpc~ 3 . If these sources produce UHECRs, 
we have P C r(10 17 - 10 21 eV)/P 7 < 0.1 for j g ~ 2.0. This 
suggests that either a small fraction of the BL Lac ob- 
jects in the distant universe are capable of accelerating 
protons to ultra-high energies or the baryon-electron ra- 
tio in these sources is significantly smaller than one. 

Though there are few astrophysical sources that can ac- 
celerate protons to energies beyond 10 20 eV, there should 
be more sources in the universe that can accelerate pro- 
tons to a lower E p ^ max . For example, accretion shocks 
in clusters of galaxies may be able to accelerate protons 
to 10 18 - 10 19 eV (e.g. Inoue et al. 2005). These ul- 
trahigh energy protons produce diffuse gamma-ray emis- 
sion as well. Therefore, the observed extragalactic dif- 
fuse gamma-ray background also constrains the UHECR 
emissivity produced by these lower E PjTnax accelerators. 
The upper limits of UHECR emissivity for different 
Ep.max are shown in Fig. 3. It shows that for lower E p>max 
accelerators, the UHECR emissivity could be higher, but 
it must still be lower than ~ 10 47 ergMpc _3 yr _1 at z = 1. 

3.4. The effect of intergalactic magnetic field on the 
cascade radiation 

When the intergalactic magnetic fields are strong, the 
synchrotron cooling of the secondary electrons produced 
by UHECRs propagating in the intergalactic space could 
dominate over the inverse- Compton cooling and as a re- 
sult, the cascade energy is reduced to some extent com- 
pared to the case without the magnetic fields. Unfor- 
tunately, very little is known about the origin, spatial 
configuration and amplitude of the intergalactic mag- 
netic fields. Only upper limit on the intergalactic mag- 
netic fields is obtained from Faraday rotation measure- 
ments, i.e. B < 10 _8 G(/ c /lMpc) _1 / 2 , where l c is the 
coherence length scale of the fields (e.g. Ryu et al. 
1998). Following Gabici & Aharonian (2005) and Kotera 
et al. (2011), the effective inverse Compton cooling 
time on the CMB and radio backgrounds can be writ- 
ten as t ei ~ 5 x 10 14 s(£ e /10 18 eV) Q ' c , with a IC =l 
if the electron energy E e < 10 18 eV and ajc=0.25 if 
10 18 eV < E e < 10 20 eV. Above 10 20 eV, t ei ~ 1.6 x 
10 15 s(-E e /10 20 eV). In comparison, the synchrotron cool- 
ing time is ( eB - 4 x 10 14 s(B /lnG^GEe/lO^eV)- 1 . 



The opposite scalings of i e7 and t e B with electron energy 
imply the existence of a cross-over energy E eyCr , above 
which electrons cool mainly via synchrotron radiation 
instead of undergoing an inverse Compton cascade, i.e. 
£ e cr = 10 18 eVS-(i for B > InG, E e , cr = lO^eVB"^ 6 
for InG > B > O.lnG, and E^ cr = 5 x lO^eV-B" 1 ^ 
for B < O.lnG (Kotera et al. 2011). Since the energy 
of secondary electrons produced by UHE protons is a 
fraction of 1/20 or 10~ 3 of the parent proton energy in 
the photopion or photopair process, electrons with ener- 
gies above E e ^ cr > 10 18 eV can be produced only by the 
photopion process of UHE protons. In Fig. 4, we show 
the fraction energy loss of UHE protons into the cascade 
EM component (i.e. excluding the synchrotron radia- 
tion) in the presence of the intergalactic magnetic fields 
(for three cases of B = 0.01, 0.1 and 1 nG ). It shows that 
this fraction drops significantly as the magnetic fields in- 
crease for protons above the photopion threshold energy. 
We further calculate the cascade energy density uj^ in 
the presence of intergalactic magnetic fields. In Fig. 5, 
we show the ratio of the cascade energy densities in the 
presence of a magnetic field (w^) and in its absence (w 7 ) 
as a function of the maximum proton energy for different 
UHECR source populations. When E p ^ max is large, up to 
~ 45% of the cascade energy is lost into the synchrotron 
radiation. 

The characteristic energy of the synchrotron radiation 
is E ltSyn ~ 6.8GeVB„ G (.E e /10 19 eV) 2 . So when the mag- 
netic fields are stronger, e.g. InG < B < lOnG, the syn- 
chrotron peak falls at 10-100 GeV. Thus the synchrotron 
component in fact also contributes to EGRB, although 
its spectrum is different from that of the cascade radi- 
ation. Therefore we conclude that the presence of in- 
tergalactic magnetic fields affects the upper limits of the 
UHECR emissivity by a factor < 2. 

3.5. The upper bound on the cosmogenic neutrino flux 

Once we know the upper bound on the UHECR emis- 
sivity, we can get the upper bound on the cosmogenic 
neutrino flux. As the energy loss distance for protons of 
energy above the photopion production threshold in in- 
teractions with CMB photons is relatively short (< 300 
Mpc), it is reasonable to assume that cosmic rays at these 
energies lose all of their energy locally. Following Engcl 
et al. (2001) and Yuksel & Kistler (2007), we use the ap- 
proximation that the fraction of the original proton en- 
ergy that is lost to neutrinos can be parametrized with 
a gradual step function ip(E p ) = 0.45/(1 + (E t /E p ) 2 ), 
where 0.45 is the asymptotic fraction of injected cosmic- 
ray energy transferred to neutrinos above 10 21 eV and 
E t ~ 2 x 10 20 eV. Thus, the total + energy flux at 
Earth can be written as a sum of the contributions by in- 
dividual UHECR sources that are generated at different 
cosmological epochs, i.e. 

^* = iF'T^* (£ ' (1 + 2))i ' ((, §l < "' 

(16) 

where E s p = 20(1 + z)E v (approximating each daugh- 
ter neutrino receiving about 1/20 of the injected proton 
energy) . 

The results are shown in Fig. 6 for Ep^min — 

10 17 eV, 

E m ax = 10 21 eV and a flat proton spectrum with 7 g = 2. 
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In contrast to Berezinsky et al. (2011) and Ahlers et al. 
(2010), we relax the requirement that the accumulated 
UHECR flux from a homogenous population of sources 
in the whole universe fit the observed energy spectrum 
of UHECRs since UHECRs above 10 19 eV are produced 
dominantly by local sources within hundreds of Mpc and 
the UHECR emissivity in the distant universe could be 
irrelevant. As a result, the upper bound on cosmogenic 
neutrino flux shown in Fig. 6 represents a true upper 
bound that is independent of the unknown density distri- 
bution of nearby sources that contribute to the observed 
UHECR flux. This upper bound is below the sensitivity 
of Auger and is marginally reachable by Icecube and the 
next-generation detector JEM-EUSO. 

The commonly-used Waxman-Bahcall bound (Wax- 
man & Bahcall 1999) for high-energy neutrinos also as- 
sumes that UHECR emissivity in the distant universe is 
connected with the flux of observed UHECRs that are 
produced by local sources. Without any independent 
constraints, the UHECR emissivity in the distant uni- 
verse could in principal be arbitrarily high. However, 
with the upper limit of the cascade radiation, the con- 
straints on the UHECR emissivity in the distant uni- 
verse becomes possible now and hence the neutrino upper 
bound becomes solid. 

4. THE IRON NUCLEI CASE 

The composition of UHECRs remains disputed. Al- 
though HiRes observations favor proton composition 
(Abbasi et al. 2010), recent observations by Pierre 
Auger Observatory (PAO) show a transition in the max- 
imum shower elongations < X max > and in their fluc- 
tuations RMS(X max ) between 5EeV and lOEeV (Abra- 
ham et al. 2010), which are interpreted as reflecting a 
transition in the composition of UHECR in this energy 
range from protons to heavier mass nuclei. However, 
one should be cautious that this claim depends on the 
poorly-understood hadronic interaction models at such 
high energies. In this section, we study the heavy nuclei 
composition case. For simplicity, a pure iron composi- 
tion above 10 19 eV is assumed in our calculation and 
the maximum energy of iron nuclei is fixed at 10 21 eV. 
As the corresponding maximum energy of one nucleon 
is only 1.8 x 10 19 eV, which is below the threshold en- 
ergy for photopion production of protons interacting with 
CMB photons even at high redshifts, we can neglect the 
photopion energy loss for these UHE iron nuclei. 

For an Fe nucleus generated at cosmological time t(z), 
it suffers from both energy loss due to Bethe-Heitler pro- 
cess and nucleon loss due to photo-disintegration during 
the propagation in the intergalactic space. Let's denote 
7jv(i) as the Lorentz factor of the nucleus and A(t) as 
the mass number of the nucleus. As the Lorentz factor 
of the nucleus is conserved during photo-disintegration, 
the energy loss is due to Bethe-Heitler cooling and adi- 
abatic expansion of the universe. The energy loss due 
to de-excitation of nuclei following photo-disintegration 
interactions is found to be always less efficient than the 
Bethe-Heitler energy loss (Aharonian & Taylor 2010), so 
it can be safely neglected. So the the evolution of 7jv(i) 
with time is given by 



dt 



where 

7N,BH(t, A) 
Z 2 r r-oo 



^A%.BH(t) 



~2^L th de 7 (J SH(e 7 )/(£ 7 )e 7 / e7/2r rfe- 



r(<0 
7T- 



(18) 



is the Bethe-Heitler energy loss rate for nucleus of charge 
Z and mass number A, j p .BH{t) is the Bethe-Heitler en- 
ergy loss rate for protons of the same Lorentz factor j^. 
The photo-disintegration results in nucleon loss of nuclei, 
so the mass number evolves with time as 



dA(t) 
dt 



= R A (t,^ N ), 



(19) 



where 



J?a(*,7jv) 



27^ 



/ dsa A (s)s I dxx~ 

Jsth Je/2j N 



2 n 1 {x) 



(20) 

is the total photo-disintegration rate for nucleus with 
Lorentz factor 7^ and mass number A, <ja{z) is the total 
photodisintcgration cross section, and £ t h is the thresh- 
old energy of the photon in the nucleus rest frame. 

It is found that the single nucleon loss is the dominant 
channel for photo-disintegration of heavy nuclei (Pugct 
et al. 1976), so as a good approximation, we here only 
consider single nucleon loss channel in the following cal- 
culation. This approximation results in an error less 
than 30% for nucleus with Lorentz factor smaller than 
2 x 10 10 (corresponding to energy < 10 21 eV for an iron 
nucleus) (Pugct et al. 1976). The cross-sections for pho- 
todisintcgration in the energy range £th < £ < 30 MeV 
with single nucleon loss is dominated by the giant dipole 
resonance (GDR), which can be approximately described 
by a Lorentzian form (Puget et al. 1976; Anchordoqui 
et al. 2007) as 



a A{e) = 73 



a 0,A£ 2 A-Q DR 



c2A2 - 
fc GDR 



(21) 



= -y N (t)H{z) + j N ,BH{t, A), 



(17) 



where Aqor arid ao,A are the width and maximum 
value of the cross section, £0 is the energy at which 
the cross section peaks. Fitted numerical values are 
a . A = 1.45A x 10~ 27 cm 2 , A GDR = 8 MeV, and 
£0 = 42.65^-°' 21 (0.925A 2 - 433 ) MeV for A > 4 (A < 4) 
(Karakula & Tkaczyk 1993). 

Combining Eq.(18) and Eq.(20), one can obtain 7jv(i) 
and A(t). Denote t A =i as the time when the mass num- 
ber of a parent nucleus drops to A = 1. Then, one can 
calculate the energy lost into the EM component by an 
Fe nucleus with an initial energy E S N during the whole 
period from the the injection time t(z) to the present 
time 

^(E' N ,t)E' N ee J^ z ] i [y N ,BH(t,A)A(t) + (56 - A) 
xip.BH{i)\m p c 2 dt + J A=1 [56%,BH{t)]m p c 2 dt, 

(22) 

where the second term on the right side is the Bethe- 
Heitler energy loss for those secondary nucleons that have 
been already disintegrated from the parent nucleus. Note 
that, for the sake of analytic calculations, we have as- 
sumed that all the secondary nucleons that have been 
disintegrated from the parent nucleus have an averaged 
7p at some time t, which is a good approximation since 



() 



the photo-disintegration process of nuclei in the major- 
ity of the relevant energy range is faster than the Bethe- 
Heitler cooling (Puget et al. 1976; Stecker & Salamon 
1999; Ave et al. 20050 In Fig.7, we show /3f m as a 
function of the energy of iron nucleus generated at dif- 
ferent redshifts. Note that the bumps in these curves 
correspond to the energies of those iron nuclei that in- 
teract with CMB photons at the GDR peak. The rise of 
f3^ m at higher energies is due to the contributions by the 
nucleons that have been disintegrated from the parent 
nucleus. For an iron nucleus generated at z > 1 with an 
energy > 5 x 10 19 eV, more than a half of the nucleus 
energy is lost into the EM cascade radiation. Compared 
with the proton case, the energy loss fraction is smaller 
at energies Ep e < 10 20 eV. 

Similar to the proton case, we can calculate the cor- 
responding energy that remains at the present epoch, 
Po.em^N^ by taking into account the redshift energy loss. 
Finally, we obtain the cascade energy density produced 
by the whole UHECR source population in the universe 



P"(E s N ,z(t))L N (z = l)- 



S(z) dt 



S(z = 1) 
( 23 ) 

where L^{z = 1) is the emissivity of UHE iron nuclei 
per energy decade at z = 1. Requiring w 7) jy ^ w ™t X : 
we obtain the upper limits on the UHECR emissivity 
at redshift z = 1 (i.e. P(z = 1)) for different source 
evolution scenarios, which are shown in Fig. 8. Note that 
for the iron composition case, P(z = 1) is defined as the 
integral of Ljy(z = 1) over the energy range from 10 19 eV 
to 10 21 eV. These upper limits are a factor of 3-5 times 
higher than that in the proton composition case due to 
smaller energy loss fractions ^ em at energies Ep e < 10 20 
eV. 



5. SUMMARY AND DISCUSSIONS 

The observed UHECR intensity above 10 19 eV reflects 
only the UHECR emissivity in the nearby universe due 
to that these extremely high-energy particles lose en- 
ergy quickly while propagating in the universe. Thus, 
there is no direct information about the UHECR emis- 
sivity in the distant universe. In this paper, we sug- 
gest that the cascade gamma-ray radiation initiated by 
the secondary particles produced by UHECRs provides 



a useful probe of the UHECR emissivity in the distant 
universe, since this cascade gamma-ray radiation can be 
observed directly. Using the new Fermi/LAT measure- 
ment of the EGBR, we obtained upper limits on the 
UHECR emissivity in the distant universe through an 
analytic treatment. Both proton composition and pure 
iron composition cases are studied. These limits are 
then compared with the gamma-ray emissivity of can- 
didate UHECR sources, i.e. GRBs and AGNs, at high 
redshifts. We find that, for a flat proton spectrum, the 
upper limit of the UHECR emissivity is a few tens times 
larger than the observed gamma-ray emissivity in GRBs, 
while it is about one order of magnitude smaller than the 
gamma-ray emissivity in BL Lac objects. We also find 
that the presence of intergalactic magnetic fields has in- 
significant effect on the derived upper limit. Further- 
more, an upper limit on the cosmogenic neutrino flux is 
obtained from the upper limit of UHECR emissivity. In 
contrast to Berezinsky et al. (2011) and Ahlers et al. 
(2011), we relax the assumption that the accumulated 
UHECR flux from a homogenous population of sources 
in the whole universe fit the observed UHECR energy 

sjgjet^um, since UHECRs above 10 19 eV originate from 

dzlocal sources within hundreds of Mpc and the UHECR 
emissivity in the distant universe could be irrelevant. 
Therefore our upper limit on the cosmogenic neutrino 
flux represents a true upper bound that is independent 
of the unknown density distribution of the nearby sources 
that contribute to the observed UHECR flux. This up- 
per bound is below the sensitivity of Auger and is only 
marginally reachable by Icecube and the future detector 
JEM-EUSO. 

Recently, Neronov & Semikoz (2011) obtained the flux 
of the EGBR using a different method from that used 
in Abdo et al. (2010), which is a factor of ~ 2 below 
the flux obtained in Abdo et al. (2010) in the 100-200 
GeV range. If this flux is true, the upper limits on the 
UHECR emissivity as well as the upper bounds on the 
cosmogenic neutrinos will go down by the same factor of 
~2. 
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E p (eV) 

Fig. 1. — The fraction energy loss of UHE protons into the electro-magnetic component for protons generated at different redshifts. The 
dashed lines and dash dotted lines represent the fractions of energy loss through the photopair channel and photopion channel respectively, 
while the solid lines represent the sum of them. 
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Fig. 2. — Upper limits on the UHECR emissivity as a function of the spectral index of the cosmic ray proton spectrum for different 
source density evolution scenarios. The upper panel shows the upper limits on the UHECR emissivity at z = 1, and the bottom panel 
shows the upper limits on the UHECR emissivity at the present time (z -C 1) if the source density evolution follows S(z) up to the present 



time. The maximum energy is fixed at E p>ri 



10 eV, while three values are assumed for the minimum energy. In the upper panel, the 



gamma-ray emissivity of candidate UHECR sources (GRBs and BL Lac objects) at z = 1 are plotted for comparison. 
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Fig. 4. — The fraction energy loss of UHE protons into the electro-magnetic component for protons generated at two different redshifts 
in the presence of intergalactic magnetic fields. Different lines correspond to different strength of the intergalactic magnetic fields. 




Fig. 5. — The ratios of the cascade energy densities in the presence and in the absence of intergalactic magnetic fields. The blue and red 
lines correspond to the ratios for intergalactic magnetic fields of O.f nG and 1 nG, respectively. Ep^min — f0 17 eV and 7 9 = 2 are used in 
the calculation. 
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FlG. 6. — Upper limits on the cosrnogenic neutrino flux (V M H-f^) , derived from the upper limits of the UHECR emissivity in the universe, 
for different source density evolution scenarios. Ep,min — lO^eV, Ep,max — 10 21 eV and 7 9 = 2 are used in the calculation. 
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Fig. 7. — The fraction energy loss of UHE iron nuclei into the electro-magnetic component for iron nuclei generated at different redshifts. 
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Fig. 8. — Upper limits on the UHECR emissivity as a function of the spectral index for a pure iron nuclei composition of UHECRs in 
the energy range from E m i n = 10 19 eV to E max = 10 21 eV. 



